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Abstract 



C , Stationary structures in a classical isotropic two-dimensional continuous Heisenberg 

ferromagnetic spin system are studied in the framework of the (2 + l)-dimensional 
Landau-Lifshitz model. It is established that in the case of S{r, t) = S{r — vt) the 
Landau-Lifshitz equation is closely related to the Ablowitz-Ladik hierarchy. This 
relation is used to obtain soliton structures, which are shown to be caused by joint ac- 
tion of nonlinearity and spatial dispersion, contrary to the well-known one-dimensional 
^1^ , solitons which exist due to competition of nonlinearity and temporal dispersion. We 

(^ ' also present elliptical quasiperiodic stationary solutions of the stationary (2 + 1)- 

0^ , dimensional Landau-Lifshitz equation. 

o 

^ ■ 1 Introduction 



As it is known, despite of the fact that magnetism is an essentially quantum effect, a wide 
range of magnetic phenomena can be successfully described in the framework of classical 
models. One of the most widely used of such models is the one by Landau and Lifshitz, 
K> ' when magnetic is considered in the continuous limit and interaction between magnetic 

C^ ■ dipoles is taken into account in terms of some effective magnetic field. The simplest case 

of the Landau-Lifshitz model is the case of the so-called isotropic continuous Heisenberg 
ferromagnetic spin system, which is governed by the equation 

dtS = g\Sx^§\, 52 = L (LI) 

Here dt = d/dt and A is the two-dimensional Laplacian, A = dxx + dyy. This equation 
attaches much attention not only from the viewpoint of its application in the physics of 
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magnetic phenomena, but also from the viewpoint of the theory of integrable nonhnear 
partial differential equations. It is known that in the (1 + l)-dimensional case, 

dtS = g[Sxd,j], (1.2) 

this equation, which has been discussed in a large number of publications (see jJEj as well 
as the books 01^151 and references therein), can be solved using the inverse scattering 
transform (1ST). Another well studied reduction of 1)1.1(1 is the static two-dimensional case, 
which may be referred to as (0 + 2)-dimensional one, 



S xAS 



0, (1.3) 



(see, e.g., [4j) and which is closely related to the elliptic sine-Gordon model. 

The subject of the present paper are the stationary structures of the isotropic two- 
dimensional classical continuous Heisenberg spin system and we look for solutions of (|1.1|) 
which are of the form 

S = S (x -Vxt,y -Vyt) (1.4) 

(the so-called Tijon-Wright ansatz with zero frequency |Sj). 

Of course, this reduction (like any other ansatz) is a necessity, if we want to proceed 
analytically, and is due to the fact that we cannot at present integrate the original (2 + 1)- 
dimensional equation. 

On the other hand, the stationary structures which are discussed below are of much 
interest for the physics of magnetism and nonlinear physics in general since they are 
realization of the so-called dynamical solitons. In physics of magnetic phenomena there 
exist two types of localized structures. First is the domain walls (kinks) which connect 
two different ground states and which cannot be destroyed without remagnetising regions 
of macroscopic sizes (that is why they are called 'topological solitons'). Another type is 
solitons, which can be viewed as bound states of a huge number of magnons. The question 
of their existence and stability, contrary to the case of the topological solitons, is not so 
trivial. In some sense, they exist due to the presence of some conserved quantities such 
as number of magnons, total momentum and energy [3j. These are the only physical 
constants of motion which do not depend on the model we use (as to an infinite number of 
conservation laws appearing, e.g., in (1 + l)-dimensional Landau-Lifshitz equation, they 
seem to be an attribute of the model and do not survive when we move to a more realistic 
one). The moving stationary structures are the simplest field configurations for which all 
of the constants are non-zero, i.e. they are the simplest of general (or non-degenerate) 
ones. Studying these structures one can explicitly see the interaction of the nonlinearity 
and dispersion which is known to be the core mechanism of creation of localized objects 
(solitons) not only in magnetic systems but in many other areas of the nonlinear physics. 

After substitution (|1.4() equation (|1.1|1 becomes 



SxAS 



+ (v,VS)=0. (1.5) 



Introducing the variables 

z = — [x-Vxt + i{y - Vyt)] , z = — [x -Vxt-i{y - Vyt)] (1.6) 
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where v = \v\, we can rewrite it as 

S X 5,,-] + X^S, + X-^S, = (1.7) 

or, equivalently, as 

S^2 + [Sj,) S+ [(^X^S, + X-^S,) xs\=0 (1.8) 

with A = exp(i7/2), where the angle 7 is defined by Vx = UC0S7, Vy = T;sin7. Equa- 
tion (|1.7|1 is known to be integrable (its zero-curvature representation (ZCR) one can find 
in the paper [S] ) , and one can tackle it by elaborating the corresponding inverse scattering 
transform. However, in the present paper we do not discuss this question. Our aim is 
to establish the relations between the model considered and the other integrable models, 
which will provide us with a wide range of physically interesting solutions. 

For our further purposes it is convenient to rewrite (|1.8)1 in the matrix form using the 
correspondence 

3 

S={Si,S2,S3)^S = J2Sa<T''^ (1.9) 

0=1 

where a" are the Pauli matrices 

i_/Ol\ 2_ f -i \ 3 f 1 



Equation (|1.8jl then can be presented as 

Sz2 + \ (tr 5,5,-) 5 + ^ [A^^, + A-^S, , 5] = 0. (1.11) 

Namely this equation is the central object of our investigation. 

We will use the following remarkable fact: equation (|1.11|) is gauge equivalent to the 
0(3, 1) nonlinear a-model [Jj in the similar way as, e.g., the (1 -|- l)-dimensional classi- 
cal continuous Heisenberg ferromagnetic spin system H1.2|) is equivalent to the nonlinear 
Schrodinger equation (see [SE]), or the Ishimori magnetic jHj - to the Davey-Stewartson 
system (see [HI). This equivalence can be briefly described in terms of the 1ST as follows: 
some combinations of the Jost functions of the linear problem associated with the 0(3, 1) 
cr-model solve equation (|1.11|) (below we shall discuss this question more comprehensively). 
The fact that model (|1.8|) . or (|1.11|) . is related to the 0(3, 1) nonlinear cr-model is a genera- 
lization of the already known result (see, e.g., j^) that in the static case equation (|1.3() is 
gauge equivalent to the elliptic sine-Gordon equation. The 0(3, 1) cr-model, as it has been 
shown in [7j, is, in its turn, closely related to the Ablowitz-Ladik hierarchy (ALH) jlUj . 
So we shall establish the direct links between the model considered and the ALH, which 
is much more well-studied than equations (|1.5j) , (|l.llj) or models [HI Ej • 

In the present paper we first derive the gauge equivalence between the Heisenberg 
equation and the ALH (Sections ElOI) and then use it to obtain soliton solutions (Section^} 
and the elliptical quasiperiodic ones (Section E)). 



Stationary Structures in 2D Continuous Heisenberg Ferromagnetic Spin System 259 

2 The Heisenberg equation and the Ablowitz— Ladik 
hierarchy 

The method used here, which may be called the 'embedding into the ALH' method, has 
been discussed in [7| lll[[T^ . Its main idea is that some equations can be, in some sense, 
'derived' from the system of differential-difference equations (DDE) belonging to the ALH, 
which means that any common solution of several equations from the ALH also solves the 
equation we are dealing with. Relatively to the problem considered this can be briefly 
outlined as follows. 

Consider the system of two equations from the ALH, 

idxqn = Pn (qn+l + Qn-l) , (2.1) 

dyQn = Pn (qn+l " Q'n-l) , (2.2) 

where 

Pn = l- qnTn, r^ = -KQn, K = ±1. (2.3) 

Equation 1)2. 1|) is the well-known discrete nonlinear Schrodinger equation (DNLSE) |13j . 
modified by the substitution qn -^ qnexp{2ix), while the next one, (|2.2|1 . is the discrete 
modified KdV equation (DMKdV) ^1]. These equations can be rewritten in terms of the 
complex variables z = x + iy, z = x — iy as 

idqn=Pnqn+i, (2.4) 

idqn=Pnqn-i, (2.5) 

where d stands for d/dz and d for d/dz. 

It is very important that these equations are compatible, since they belong to the same 
hierarchy, and the constants of motion that play the role of the Hamiltonians for the flows 
(|2.H) . (|2.2)) are in involution. Hence, we can consider them simultaneously, as one system 
of two equations. It has been shown in 7 , and one can easily verify this fact by simple 
calculations, that, for any fixed n, each solution of system (|2.1|) . (|2.2)1 also solves the field 
equations of the 0(3, 1) nonlinear ir-model, 

dBq+^^^^^^^ + {l-qr)q = Q, r = -nq (2.6) 

1 — qr 

and that the quantities pn satisfy the 2D Toda lattice equations 

dd lupn = Pn-l - '^Pn + Pn+1 (2.7) 

(see ^Jl). Namely this we bear in mind when say that the 0(3,1) nonlinear cr-model and 
the 2D Toda lattice can be 'embedded' into the ALH. 

The situation with the Heisenberg spin system is somewhat more difficult. Solution for 
equation Hl.ll() cannot be constructed by means of g„'s and r„'s only. To do that we have 
to consider the ZCR for the ALH and to analyze the corresponding linear problems. 



Un = 




Tn \ 
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The integrable DDEs (|2.4j) . (|2.5j) . as well as all equations of the ALH, can be presented 
as the compatibility condition for the linear system 

*„+l(A) = Un{\)^n{\). (2.8) 

d^nW = K(A)^„(A), (2.9) 

d^n[\) = ^„(A)^„(A), (2.10) 



(2.11) 



Vn = -i\ : V "'' 1, Wn=i( ^' ^""'''" ^'^^ ). (2.12) 

y X '^Qn A ^ - rn-iqn J \ Mn-l J 

One can easily see that equations (|2.8j) . (|2.9|) and (|2.8|) . (|2.1()|) are compatible only if 
matrices C/, V and T^ satisfy the so-called zero-curvature equations 

dUn = Vn+lUn " t/nK (2.13) 

and 

BUn = Wn+lUn " C/„T^„ (2.14) 

which are equivalent to (|2.4|) and (|2.5|) correspondingly. 

Namely the solutions of the linear problems H2.8() - ()2.10() . ^„,'s, are the key objects of 
our consideration and the main result of this paper can be formulated as follows: for any n, 
matrices 

Sn = ^-^a^n (2.15) 

constructed of solutions of the linear problems of the ALH solve the matrix Landau- 
Lifshitz equation (|1.1H) . 

To derive this result consider matrices o"^, o = 1,2,3 defined by 

< = < (A; z, z) = ^-\\- z, z)a'^^„(A; z, z), (2.16) 

where a'^ is a Pauli matrix (|1.10j) . and ^„, recall, is a matrix solution of system (|2.8j) . 
(tmH) (in this notation 5„ = al). It follows from (OTItI) and ^I^, (tmH) that 

dal = ^-1 K, K] ^n, 9< = ^-1 K, VF„] ^„. (2.17) 

Using expressions (|2.12|) for Vn and Wn one can find the derivatives of the matrices Sn in 
terms of the matrices S^ given by 5^ = (1/2) (cr^ it ia'^ as follows: 

(iA/2) dSn = r„_i5+ - qnS-, (2.18) 

(i/2A) dSn = -rnS+ + Qn-iS-. (2.19) 
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These relations together with analogous expressions for the derivatives da^^ Ba^ and 
formulae H2.4|) . H2.5|l . after straightforward calculations, omitted here, lead us to 

-ddSn = {qn-irn-i + qnTn) Sn + (A"V„ - Ar„,_i) 5+ + (A"^g„_i - Xqn) S~ . (2.20) 
Noting that 

tr dSn BSn = -4 {qn~ir„^i + qnrn) (2.21) 

and 

[X^dSn + X-''dSn , Sn] = 4i (Ar„_i - A-V„) S+ + 4z (Ag„ - A-^(7„_i) 5" (2.22) 

(both of these formulae follow from ()2.18|) . ()2.19|) ) we obtain that for every n the matrix Sn 
solves the equation 

ddSn + \ (tr dSndSn) Sn + J. [X^dSn + X'^BSn , Sn] = (2.23) 



which is the main equation of our study, (see (jLllj) '). This key result of the present paper 
can be reformulated in terms of the vector Sn, which corresponds to the matrix 5„ and 
which can be presented as 

(2.24) 



Sn — {Snl, Sn2, Sns) , 






where 






*'"'*•" 


Sn3 - 


*^i^^)^^r)+^^^^^^^^^) 


'"' '" *<"'*<?'-*<"'*<?"• 


^^^^^^i^^^-^i^^)^^' 



(2.25) 

(here ^n are the elements of the matrix V&n), as follows: for each n the vector 5„ defined 
by ((2211), (p^ solves equation (fTsll . 

Thus we have established the links between equation (jLllj) . or ()1.8p . describing sta- 
tionary moving structures in the (2 + l)-dimensional classical continuous Heisenberg spin 
system and the ALH. Some more detailed analysis of the gauge equivalence between these 
models one can find in the next section. However, in this work we are going to focus our 
attention on 'practical' aspects of this relation, so a reader can consider it as an 'empirical' 
fact which can be straightforwardly, and rather easily, verified by the calculations outlined 
above. 

As was mentioned earlier, model ()1.7() is known to be integrable and its zero-curvature 
representation has already been written out. But, to our knowledge, the corresponding 
1ST has not been elaborated yet, while the ALH is one of the best-studied nonlinear in- 
tegrable models. Besides, the Heisenberg equation is a vector problem, which somehow 
complicates inverse scattering analysis, while the ALH is a scalar one. So, to our opin- 
ion, the 'embedding into the ALH' approach is rather promising and in what follows we 
demonstrate its usefulness by constructing the soliton and quasiperiodic solutions for the 
equations considered using the already known solutions for the ALH. 

The magnetic energy density, W, 



W 



S 



I (VS, VS) (2.26) 



262 G M Pritula and V E Vekslerchik 

of the field configurations obtained by the embedding into the ALH method can be ex- 
pressed in terms of the Qn and r„'s: 



W 



^n 



V^ 



-7- [qnTn + Qn-irn-l) ■ (2.27) 

45 



It can be shown that from the viewpoint of application of solutions of the ALH equations 
to the description of the vector field S one has restrict himself with the case of k = 1 , 

rn = -Qn, Pn = l + kn|^ (2.28) 

when the components of the vector 5„ (|2.24|) are real (in the opposite case, k = — 1, the 
components of S„ are complex) and the magnetic energy ()2.27|) is positive. 

In the next section we will consider the relation between equation (|1.1H) and the ALH 
in the framework of the 1ST. 

3 Gauge equivalence and zero curvature representation 

In the previous section we considered the relation between the ALH and the Landau- 
Lifshitz equation in terms of solutions: we demonstrated how to use solutions of the ALH 
to obtain ones for the Landau-Lifshitz equation. Now we are going to discuss this question 
in somewhat more general way. Both the ALH and the Landau-Lifshitz equations are 
integrable models and it is interesting to describe this correspondence in the language of 
the 1ST and to derive links between the auxiliary linear problems which are used to present 
the integrable models in the zero curvature form (namely this is usually understood when 
one uses the words 'gauge equivalence'). 

Let us consider again the auxiliary linear problems of the ALH mentioned in Section [2l 
To our current purposes we do not need the discrete problem 1)2. 8() and will be dealing 
with the continuous ones (|2.9|) . ()2.1U|) . So, we omit now the index n and rewrite (|2.9() . 
(fTT71|) . (tTT^ as 

9^(0 = ^(C)^(C), 9^1/(0 = wiO^iO, (3.1) 

where 

\C 'qi C - roqi J 
and 

.(C)..(^=--^-) ,3.3) 

(we have replaced qn, rn with qi, ri and qn-i, '"n-i with qo, vq). In what follows we 
denote the spectral parameter by C and use A for its particular value appearing in the 
definition (|2.1Hj) of the matrix S, 

S = S{z,z) = ^"^(A; z, z)ct^*(A; z, z) (3.4) 
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The compatibility (zero-curvature) condition for the system (|l-{.lj) 

dV{C)-dW{C) + [V{C),W{C)]=0 (3.5) 

leads to the following system of four partial differential equations (PDE) for four unknown 
functions qi, ri, qo, tq: 

idqo=poqi, (3.6) 

idri = -piro, (3.7) 

idqi = piqo, (3.8) 

idro = -pori. (3.9) 

This system is in some sense intermediate between the DDEs (|2.4j) . (|2.5|) and the PDE 
1)2. 6|) : both of them can be 'reconstructed' from H3.6|) " H3.9|) (we will return to this question 
below). And namely system (|3.6|) - (|3.9j) is, strictly speaking, gauge equivalent to the spin 
field equation we are dealing with. 

Now we will derive the ZCR for the stationary (2 + l)-dimensional Landau-Lifshitz 
equation from (|3.H) using the gauge transformation by means of the matrix ^(A). Intro- 
ducing the matrix function $(C) 

$(C) = ^-Ha)^(C) (3.10) 

one can obtain from (|3.1() that it satisfies the following equations 

a$(C) = Vl{CMC), d^c) = Wl{CMC), (3.11) 

where 

VUO = ^-'(A) [F(C) - V{X)] ^(A), (3.12) 

WLiO = ^-HX) [W{C) - W{X)] M/(A). (3.13) 

Noting that 

5. = ^-\X) [a', V{X)] ^(A) = 2iX-'^-\X) ( ° '^^ ) ^(A), (3.14) 

S, = ^-HX) [<t',W{X)] ^{X) = 2iX^-\X) ( _°^ Q^ ) ^{X) (3.15) 

one can present Vl and Wl as 

Vl(C) = ^ (C-' - A-^) iS-l) + l {XC' - 1) SS,, (3.16) 

Wl{0 = ^ (C - A^) iS + l) + \ iCX-' - 1) SS,. (3.17) 
Using the zero-curvature conditions for equations (|3.11|) . 

BVl - dWL + [Vl, Wl] = 0, (3.18) 
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and calculating the left-hand-side part of this equation 

dVL-dWL + [VL,WL] = 

1 (^ C\ \ aa , 1 /o o . a a \ . ;\2q , ■\-2, 



one can conclude that the matrix S must solve 

S5, + l(s,s, + 5.5.) + a^s, + a-s.- = o. (3.20) 

Noting that the anticommutator of traceless 2x2 matrices is proportional to the unit one 
and that the anticommutator of S and Sz or S2 is zero (which follows from the fact that 
det 5 = 1, which is another form of the equality S"^ = 1) one can present this equation in 
the form (jLllj) . Thus the linear problems p.llf) together with definitions ()3.16|) and ()3.17|1 
can be viewed as the ZCR for the main equation of the present paper. 

After we have derived the ZCR for Hl.ll[) we would like to make a few remarks on 
the application of the inverse scattering technique to non-evolutionary equations as ours. 
The 1ST has been originally developed for the Cauchy problems. However, since then 
much efforts has been made to adjust this method for various boundary value problems. 
This is a rather difficult task since the latter seem to be more difficult than the former 
ones. Among successes in this field one should mention results related to the hyperbolic 
systems such as, e.g, the sine-Gordon equation, the principal chiral field equations etc. For 
these models the initial value - boundary value problems has been shown to be well stated 
problems, the existence and uniqueness of the solution has been established and IST-based 
algorithms to solve, say, the Goursat type boundary problem have been elaborated. One 
can find discussion of some recent results on boundary problems for the (l + l)-dimensional 
systems on semi- infinite and finite interval, for example in |15j . 

As to the elliptical systems, similar to the one discussed here, which do not possess 
characteristics, it is also possible to apply the 1ST for solving some boundary value prob- 
lems. Usually it is achieved by breaking the symmetry between the coordinates (which 
is, of course, not very natural for this kinds of equations), selecting one of them (say, y), 
considering the problem on a half-plane {y > 0) or on a finite domain (in this case the 
part of the boundary data plays role of the Cauchy conditions) and performing analysis 
(i.e. solving the direct and inverse spectral problems) for the auxiliary linear equation cor- 
responding to the complementary coordinate (say, ^x = U^). One can find examples of 
such approach in 16 and references therein. However, in this paper we do not discuss the 
mathematically rigorous formulation of the problem related to (|1.8|) . We consider here the 
1ST as a method to generate some classes of particular solutions and restrict ourselves to 
the ones most interesting from the physical viewpoint, solitons and quasiperiodic solutions. 

Above we have mapped the V-W pair for system (|H.6j) - (|:i9|) into the Vl-Wl pair for 
equation ()1.11|) by means of the gauge transformation 1)3. 10(1 , 

VUO = -^-HX)d^{X) + *-i(A)y(C)^(A), (3.21) 

WLiC) = -^-\X)d^iX) + ^-HA)F(C)^(A). (3.22) 

Now we are going to derive the inverse transform: from (|3.11() to 1)3. If) (i.e. from the 
V-W pair (l3T6|) . (ITTTll to (|S21), ^M)- The fist step is to diagonahze a solution of the 
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Landau-Lifshitz equation, i.e., to calculate, for given S, the matrix ^ defined by 

S = ^-^a^^. (3.23) 

It is obvious that the S — > ^ correspondence is not one-to-one. For any ^ satisfying (|3.23() 
the matrix D^ with an arbitrary diagonal matrix D will also solve (|3.23|) . The main point 
of the Landau-Lifshitz equation — > ALH transform is to use this arbitrariness to present 
the matrices d"^ ■ ^~\ 5^ • ^"^ in (P0|) . (PU?|) form with C = A. 

5^.^-1 = y(A) and 9^ • 'I'"^ = l^(A). (3.24) 

This step needs some calculations which are presented in the Appendix. Performing then 
the gauge transform with the found matrix ^, one can obtain that the transformed Vl, 
Wl matrices 

y(() = 9^.^-1 + ^14(0^-1, (3.25) 

W{C) = B-^ ■ ^-^ + ^Vl(C)^~^ (3.26) 

are exactly of the form (|3.2)) . (|3.3)) which means that the functions qq, tq, qi, ri (which 
are defined now in the terms of the matrix ^ (i.e. in the terms of the matrix S) solve the 
system (|HS|-(|Sin|- 

System (|3.6() - (|3.9() that can be rewritten as the DDEs from the ALH. Indeed, starting 
from the quantities qo, tq, qi, ri one can define the quantities q2, ^2 

g2 = iz , r2 = -I- (3.27) 

1 - qiTi 1 - qiTi 

and demonstrate that they satisfy the following identities: 

idq2 = (1 - q2r2) qi, -idr2 = (1 - g2?'2) n- (3.28) 

Analogously, the quantities q^i, r_i, 

" ^ ^^° r_, = -z-^^ (3.29) 



1 - qoro 1 - qoro 

satisfy 

idq-.i = {I- q-ir-i)qo, -z9r_i = (1 - g_ir_i) tq. (3.30) 

This procedure can be repeated in both directions 

••• ^ {Q-i,r-i) ^ {qo,ro,qi,n) -^ (92,^2) ^ ■•• . (3.31) 

This gives an infinite sequence of q^s, r„'s which solve 

idqn=Pnqn+i, (3.32) 

-idrn = PnTn-l (3.33) 
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and 

idqn=Pnqn-i, (3.34) 

-idrn = PnTn+l, (3.35) 

i.e. the Ablowitz-Ladik DDEs. 

To conclude this section we want to discuss the fohowing question. If we start with 
the ALH, which is a system of DDEs, then the relation between the discrete equations 
(ALH) and the partial differential Landau-Lifshitz equation is rather obvious: our PDE is 
a differential consequence of the DDEs. But if we start with the Landau-Lifshitz equation, 
then what role do the DDEs from the ALH play in the theory of our PDE? In simpler 
words, what does the subscript n mean in terms of our PDE? The answer is as follows. 
We have an example of the situation studied by Levi, Benguria ^Z[E|, Shabat, Yamilov 
J19j and others: discrete integrable equations (the equations from the ALH in our case) 
describe sequences of the Backlund transformations for some PDEs (the stationary (2 + 1)- 
dimensional Landau-Lifshitz equation in our case). Indeed, if we have a solution of our 
equation, Si, we can derive from it the matrix ^i which solves the linear problems for the 
DNLSE and DMKdV, and hence the quantities qo, tq, qi, ri which solve (|,'-{.6|) - (|l'{.9|) . Then 

we can construct the new ^-matrix ^2=1 x _i I ^i> and the new spin field So which 

\qi \ ' J 

corresponds to the matrix ^2^ (J'^^2. This vector field will also solve the Landau-Lifshitz 

equation. This procedure can be repeated infinitely 

• • • ^ 5„ ^ *„, g„, r„ -^ ^„+i = C/„^„, -^ Sn+i -^ ■■■ . (3.36) 

Moreover, it can be performed in other direction 

>Sn^^n, Qn, r„ ^ *„-l = C/^^i^n ^ Sn-1 ^ ■ • • . (3.37) 

Thus we can obtain an infinite number of solutions ( . . . , ^-i, Sq, Si = S, S2, ■ ■ ■) from 
one solution S and relations between the vectors Sn with different values of the index n 
(Backlund relations) can be described by the equations which are analogous to (and can 
be derived from) the DDEs from the ALH. 

4 Soliton structures 

The discrete nonlinear Schrodinger equation (|2.1() under the condition (|2.28|) has been 
already solved in the pioneering work by Ablowitz and Ladik J^]. As to the solutions 
of (|2.2|) . or system ()2.4|) . (|2.5() . they can be obtained by minor modifications of the ones 
for (|2.1j) . which is again a manifestation of the fact that all of them belong to the same 
hierarchy. We will not repeat here the derivation of the 1ST (one can find the technical 
details in j^ or, say, in the book [13!) and write down only some final formulae that will 
be used below. 

The A^-soliton solution of equations 1)2. 4|) . (|2.5|) can be presented as follows: 

N N 

Uz,z)=Y,Cnj{z,mnj{z,z), ^nj{z, z) = ^^ [MnHz, z)^^^' ■ (^-l) 

j=l k=l 
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The constants A^'s are the eigenvalues of the corresponding scattering problem (|2.8j) (to 
be more precise, the discrete spectrum of the scattering problem (|2.8|1 consists of A^ pairs 
of the eigenvalues (A^, —Xk))- The functions C„fc(z, z) are given by 

Cnkiz, z) = ClXl" exp {i0fc(z, z)} , (4.2) 

where C^'s are arbitrary constants, 

^k(.z,z) = X-^z + Xlz, (4.3) 

while the matrix M is given by 

M = I + AAnAAn. (4.4) 

Here / is the N x N unit matrix, the overbar stands for the complex conjugation, 

A = diag (Ai,...,Aiv) (4.5) 

and An is the N x N matrix with the elements 

Ali.'Hz,z) = ^I^ (4.6) 

Solution for system (|2.8() - (|2.1U() can be presented in the pure soliton case as 



exp(iA^2:) 

exp(-iA^^2;) 



^n{z,z)=FniX;Z,z)[ ^\ ^ 2.^ 1' (4-7) 



where Fn{X) is the matrix of the following structure: 

1 / /n(A) gn{X) \ 

Fn{X) = y-- . _ . (4.8) 

fnioo) \ -g„(l/A) ^(1/A) J 

Here 

fn{X) = l-X^Y. '^'_;^212 ^ 5nW=AEx?^- (^-S) 

The above formulae contain all we need to construct solutions for the Landau-Lifshitz 
equation (|1.8j) . or (ll.llj) . The vertical component of the vector Sn (see (|2.25j) ) can be 
presented, using H4.7() - ()4.9() and the fact that in our case A = exp(i7/2) = 1/A (see 
remark after (|1.8|) ). as 

5.3 = 1 - ^^ (4.10) 

while the horizontal components can be written as 

Snl = Sn± COS (pn, Sn2 = Sn± sinc^„, (4.11) 



268 G M Pritula and V E Vekslerchik 

where 



Sn± = ^l-Sl-„ ipn = aTgfn{X)-aVggn{X)+X-^Z + XH. (4.12) 

The magnetic energy density H2.26|) in this case can be presented, using H2.27|) and the 
identity Pn = det F„_|_i/ det F„ as 

W|S„l = ^|^^ + %f^-2|. (4.13) 

% l/n+l(oo) /„(oo) J 

Noting that, for a fixed value of the index n, the dependence on n can be taken into 
account by the redefinition of the constants C^, we may chose n = and write the final 
formulae as follows: 



^ 2\go{e '^/^-z,z 



|2 



S,iz, z) = l- ^^ L_^ (4.14) 



and 



Si{z,z) = S±{z,z) cos ip{z, z), (4-15) 

S2{z,z) = S±{z,z) sin (p{z,z), (4-16) 

where 

S^{z,z) = ^l-Sliz,z), (4.17) 

V5(z,z) = arg/o (e''^/^;z,z)-arggo (e'^/^; z,z) + e-'''z + e'^z (4.18) 



while the distribution of the magnetic energy of the field configuration given by (|4.14jl - 
H4.16() can be written as 

>^ = ^lTr^ + V/°"'''-? -4- (4.19) 

45 l/i(oo,z,z) fo{oo,z,z) J 

These formulae describe the A^-soliton solutions of the stationary (2 + l)-dimensional 
Landau-Lifshitz equation. 

To make clear what kind of solutions we have obtained from the ALH-solitons let us 
consider in a more detailed way the simplest of the above solutions, namely the one-soliton 
ones. In this case the quantity An{z,z) = An {z,z) (see (|4.6p ) using the designation 

A? = exp(-2(5 + i7i) (4.20) 

can be rewritten as 

An{z,z) = exp{xn + iipn} , (4.21) 

where 

Xn = 2 sinh 26 (sin 71 • Re z — cos 71 • Im z) — 2n5 + x* ) (4.22) 

ijjn = 2 cosh 26 (cos 71 • Re z + sin 71 • Im z) + 77,71 + i/j^, (4.23) 
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and X*) V'* 8-1^6 some constants. Setting n = 0, returning to the real coordinates x, y and t, 
and introducing the vectors 

t;sinh2(5 / sin 71 \ - t; cosh 2(5 / cos7i \ /, „.x 

k± = , h = • > (4-24) 

2g \ —cos 71 J " 2g \ sm7i J 

one can rewrite these formulae as 

A{x, y, t) = exp {x{x, y, t) + ii;{x, y, t)} (4.25) 

with 

x{x,y,t) = ik±,f—vt] + const, (4.26) 

il){x,y,t) = (k\\,f—vt] + const (4.27) 



(here braces stand for the usual scalar product: for k = col{kx,ky), {k,r) = kxX + kyy), 
from which one can derive the following expressions for the components of the vector S. 
The vertical components, S3, can be written as 

2sinh2 2(5 1 

cosh 2(5 — cos 2r cosh 2x + cosh 25 

where F = (71 — 7)/2, or, equivalently, as 

53 = cos0 (4.29) 

with the angle 9 being given by 

-1/2 



tan - 
2 



sinF \ / cosF 

cosh X ] +1 TT ^^'^ X 



^ sinh 6 J \ cosh 6 

The horizontal components 6*1,2 can be presented as 



(4.30) 



f)=s.ne('^^'^ ] (4.31) 

S2 J \ smip ' 



with 



ip{x,y,t) = ip{x,y,t) + (ko — k^\,f — vtj + const. (4.32) 

Here, the vector ko, 

fco = f f^°^^), (4.33) 

2g \ sin7 J 

is parallel to the velocity vector, /cq = v/2g, and is related to A;_l, k\\ by fcg = fen — /c^. The 
function (p is given by 

(p = arctan (tanh 5 cot F tanh x) ■ (4.34) 
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The magnetic energy of this field configuration can be written as 

^^^^^^J_+cosh2^^osh2x^ (4.35) 

(cosh 2x + cosh 25) 

where k^ stands for |A;_|_|. 

It can be shown that W (j4.35|) is a second-order polynomial in 53: 

W = -gk\Sl + (v, fc||) 53 + Wo, (4.36) 

where Wo is some constant. The linear energy density, Wim, 



Wnn=/ dxxW(x,y,i), x^=\-^,r\ (4.37) 




can be easily shown to be 

Win = 4gfc^. (4.38) 

To simplify the following analysis let us consider the case when the velocity vector is 
directed along the x-axis (7 = 0, A = 1). This does not lead to loss of generality because 
solutions corresponding the arbitrary vector v = {v cos 7, v sin 7) can be obtained from the 
ones presented below by the substitution x — > XC0S7 + ysin7, y — > ycosj — 3;sin7. It 
can be easily seen that formulae ()4.28|) - H4.32() in the limiting cases 71 = tt/2 and 71 = 
describe essentially different field structures. In the case 71 = 7r/2 (T = n/A) 

X = x{^ ~ ^*) = k±{x — x^ — vt) (4.39) 

with an arbitrary constant x*, and both S^ and W depend on x — vt only, 

e = e{x- vt) = 2 arctan ( si^25_ \ ^ ^^^^^^ 

\[coA,\x + S) + co^\i\x-5)\' ) 
W = W{x-vt)=gkls].n^e{x-vt) (4.41) 

while if can be written as 

ip{x, y, t) = (p{x — vt) + ko{x — vt) + k\\y + const, (4.42) 

where 

0{x) = arctan (tanhJ tanh A;_l(x — x*)) . (4.43) 

So, this solution describes a localized structure, moving in the x-direction, which is phase 
modulated in the transversal direction (y-direction), and it may be termed 'quasi one- 
dimensional soliton' (see Fig. 1). 

The soliton obtained above is essentially two-dimensional structure and despite ap- 
parent similarity it cannot be reduced to its one-dimensional analogue. Indeed, in the 
one-dimensional case soliton solutions of equation H1.2() possess the following form: 

X — vt\ ^ ^ ( X — vt\ , . . .^ 
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"quasi 1 D" 



"channel' 




Figure 1. Limiting cases of the one-soliton solutions (schematically) corresponding to F = 7r/4 
and r = 0. 

and one can say that such solitons exist due to the temporal phase modulation of the whole 
medium, which manifests itself in the fact that L{Q) oc il~^'^, i.e., soliton vanishes with O 
going to zero. In other words, these one-dimensional soliton structures do not exist in 
absence of the phase modulation Qt. In our case, existence of solitons is due to the spatial 
phase modulation (in y-direction), which manifests itself in the fact that the magnetic 
energy density is proportional to k?,. In other words, solitons we have obtained differ from 
their one-dimensional analogues in the physical mechanism lying in their background: 
they are caused by the competition between the spatial dispersion and nonlinearity while 
in the one-dimensional case solitons are caused by the competition between the temporal 
dispersion and nonlinearity. 

In the opposite case, 71 = (F = 0), 



9 = 6{y) = 2 arctan 



cosh (5 



sinhA;_L(y -y*)' 



i.e. Ss, which can be written as 



^3 = S3{y) 



sinh^ k± [y — y^,) — cosh^ (^i/2) 
sinh k± {y — y^) + cosh {61 /2) 



(4.45) 



(4.46) 



depends on y only (and, what is essential, does not depend on time), while the horizontal 
components are rotating with constant frequency: 



(p = ip{x — Vt) = Kx{x — Vt), 



(4.47) 



where 



k\\ — ko 



(4.48) 



and ka, remind, stands for \ka\- 

This solution describes the spin wave localized in the 2/c^ -neighborhood of the line 
y = y* (this field distribution, which is depicted schematically in the Fig. 1, may be termed 
'channel'). The magnetic energy of the 'channel' field configuration does not depend on 
time, W = W(y), hence it can be considered as almost static, in the sense that we have 
no energy transport in this case. Similar structures have been found by A S Kovalev j2Uj . 
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Figure 2. One-soliton solution for A = 1 and Ai — 0.8exp(i7r/12) 




Figure 3. Two-soliton solution for A = 1 and Ai = A2 = 0.8exp(i7r/12) 

The character of the soliton field structure in the general case < 71 < tt/2 can be 
seen from the Fig. 2. The many-soliton structures in the general case can be viewed as 
consisting of several intersecting solitons. One can find typical two-soliton spin distribution 
in the Fig. 3. 

5 Quasiperiodic structures 



The ALH in the quasiperiodic case is less studied than in the soliton one. Several authors 
have discussed the quasiperiodic solutions (QPS) for the discrete nonlinear Schrodinger and 
the discrete modified Korteveg-de Vries equations (see, e.g., [221 Ej), but their results are 
not enough to construct the corresponding solutions for the Heisenberg equation using the 
'embedding into the ALH' method. What we need and what is absent in the papers |22 | l21 j 
is a solution of the auxiliary system ()2.8() - ()2.1U|) which in the quasiperiodic case is known 
as the Baker-Akhiezer function. Later this question was solved in !23] (see also |24j). 
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However, these results, describing general finite-genus solutions, are rather cumbersome, 
so here we restrict ourselves only to the elliptic solutions, which are the simplest QPS. 
The elliptic solutions for system ()2.4() . H2.5() possess the following structure: 

Hn — (/*c, e Q /A \ ' 'n — '*fc C q /a \ ' W'-'-J 

where ??i is one of the elliptic theta- functions (see, e.g. HSI); the phase (p is some linear 
function of the coordinates z and z (it will be specified below), 

C„ = y + ^ + n/3 + const (5.2) 

and the constants q^:, r* are related by 

^f(/3) 
^*'^* i9i(a + /3)??i(a-/?)' 

The quantities p„'s can be presented as 



(5.3) 



^Ija) ^,(C„ + /3)^,(C„-/?) 

^" T?i(a + /3)^i(a-/3) t?2(^^) i-J 

which follows from expressions (|5.1|) and the identity 

t??(x)T9i(y + 2;)T?i(y - ^) + i??(y)T?i(z + x)'&i{z - x) 

+ dl{z)Mx + y)Mx-y) = ^- (5-5) 

This identity is the Fay's formulae j26j for the elliptic functions. It can be used to calculate 
the derivatives of the ??-functions. Differentiating (|5.5|) with respect to z and putting z = y 
one can obtain for the logarithmic derivative ip, 

^(C) = ^ln^i(C), (5.6) 

the relation: 

Using the latter one can obtain that functions 1)5. 1|) satisfy equations (|2.4I) . ()2.5() provided 
the scales L and L are chosen as 

L = ,.;(0)ili^.-, Z = -i#,(0)il<^., (5.8) 

T9i(a)T?i(/3) T9i(a)T?i(/3) 

and the phase (/> is given by 

cj) = -i [iPiP) - V'(/3 + a)]y + i [V(/3) - HP - a)] 4 + const. (5.9) 
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The Baker- Akhiezer function of our problem (i.e. the quasiper iodic solution for sys- 
tem 1)2. 9p ). 1)2. lUj) can be written as a matrix with the elements 

*(i^)(A) = Al^^'l^^^^^eM^^l), (5.10) 

^rHA =-^i/"?A Qn-i '^.^.'":'^ ^ -^exp^$l , 5.11 

Tfi2)/xN . nx e{a + (3)e{r] + p)e{C„-r]-a- 13) ,.^, , ^^, 

4^2) (A) = ^2/.^^^|^^exp(.cl>2). (5.13) 

Here Ai^2 are arbitrary constants which are of no importance for our further consideration. 
The phases $1^2 are the linear functions of the coordinates, 

tz tz 

^i = [tPiv + « + /?)- V(a)] J + mv + P)+ V'(a) - V'(« - /3) - ^(Z?)] j, (5.14) 

tz tz 

$2 = [-V'(7?) + ij{a + /?) - V(a) - ^(/3)] - + [-ij{7i + a) + V(a)] -^, (5.15) 

the quantities ^1,2 are given by 

_ e{a)eirj + a) _ ^ g(a)g(r? + a + /3) 

^^ ^(a - /3)0(r/ + a + /3) ' ^^ 0(a + /3)e(r/ + a) ^ ' ' 



and 77 can be determined as a solution of the equation 

eA^ (5.17) 



Ma-P) Mv)Mv + a + (3) _ ,2 



i?i(a + /3) Mri + a)Mv + P) 

(in the framework of the general theory, rj can be considered as the point of the Riemann 
surface that corresponds to the point eX^ of the complex plane). 

These formulae (we do not present here the corresponding derivation procedure) can be 
verified straightforwardly using (|5.7() and (|5.5|) . They provide all necessary to construct 
the elliptic solutions for the Heisenberg equation (|1.11|) . Using (|2.25|) . (|4.29|) and (|4.31j) . 
and omitting the n-dependence one can obtain 

'''' 2 ^,^^ + a)Mv + a + f3) MC + V)MC-V-P) ' 

and 

exp {2i(p} = — — — exp {2i<i>} , 5.19 

^i(C -r]- /3)t?i(C -rj-a- (3) 

where 

$ = $i-$2 + </'. (5.20) 
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The last two formulae can be rewritten as 

99 = (p{C) - (p'{0)C - (K||,r-- vt) , (5.21) 

^^^^ = ^ ^" MC-,-P)MC-v-a-py ^'-''^ 

Here the vector ku is given by 

K|| ( ^°^^i 1 (5.23) 

" ' sm7i ' 



with 



2 v^ ^f{a)^f{P)^l{2r] + a + p) 



and the angle 71 is defined by 

£2 ^^(a + /5) ^ ^ 

The magnetic energy density W ()2.26|) of the above field configuration is, as in the 
one-soliton case, a second-order polynomial in 53: 

W = -g^Sl + {v, K||) 53 + Wo, (5.26) 

where Wq is some constant. The last formula again illustrates the importance of the 
transversal modulation (space dispersion) for the existence of our nonlinear structures. 

It should be noted that to ensure reality of all physical quantities, such as 5,, W one 
has to impose some restrictions on the parameters a, /? and r] (or e) which appear in 
the above expressions. We cannot at present formulate these restrictions in their general 
form, but will show below how these parameters should be chosen in some particular case, 
which is a generalization of the pure soliton one, in the sense that the one-soliton solutions 
obtained in Section |3] are some limiting cases ot the elliptical ones discussed below. 

Thus, in what follows we restrict ourselves with the case of 

a = ^, (5.27) 

where r is the complex half-period of the '!?i-function (see |13])- ^^ ^^^ be shown that in 
this case both the components of the vector S and the energy W will be real if we choose 

C = ^-iC, rj = fi-^, P = 2ip, (5.28) 

where hats indicate that correspondent quantities are real. In what follows we use together 
with the theta-functions also the Jacobian elliptical functions sn, en and dn, 

.,„ = .„(„,.) = *^^l(!^ (5.2.) 

W2(0) 't}o{u/2K) 
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with 



k 



^i(O)' 



K = ^^^(0) 



(5.30) 



(the definition of i92,3 and analogous formulae for cnu and dnu one can find, e.g., in \25\). 
The 'coordinate' C can be written now as 



c 



where 



k± 



— ik±,f—vt 



V vr sn4iK/3 f sin 71 
- cos 71 



2g 2K i 
The expressions for and ip (|5.18j) and (|5.21|1 . (|5.22|) can be presented as 



su{2Kf] + 2iKi3] 



dn 2Kfi + 2iKC 



cn(2Kfi + 2iKC 



1 " 
tan - 



and 

ip= [kq - K\\,r- vt) + arg -i?2(^? - iO'&siv - iC) 
The vector kq is given by 

Re ■00 {v + iP 



2 

Ko = - 

IT 



where 



Mz) = Hz) - K 



fci 



en 2Kz dn 2Kz 



sn2Kz 
while the vector km is given by (|5.23|) with 

V sn 4:iKf3 



V'(^) - ^^^(0) 



Mz)Mz) 
Mz)Mz) 



I — 



■ sin 2r 



g sn AKfj 

and r = (71 — 7)/2 is related to /3 and 7) by 

^2ir ^ sn2K{T] + 0) 
su2K{f}-iJ3)' 



(5.31) 



(5.32) 



(5.33) 



(5.34) 



(5.35) 



(5.36) 



(5.37) 



(5.38) 



It is straightforward to show that the limiting case of the elliptic quasiperiodic solutions 
presented above is solitons obtained in Section 01 Indeed, with the parameter k (|5.3()jl 
going to zero (which corresponds to r — > ioo), the elliptic functions sn, en and dn become 
sin, cos and 1 correspondingly. Noting that K[k = 0) = 7r/2 and identifying vr/? with 5 
(which implies vr^ — > x) one can transform H5.33() and H5.34() to formulae ()4.3U|) and (|4.32|) 
describing solutions of the Landau-Lifshitz equation in the one-soliton case. 
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6 Conclusion 

To conclude, we want to summarize the main results and to outline some perspectives of the 
studies discussed in this paper. Prom the mathematical point of view, our main result is the 
established relation between the Landau-Lifshitz equation (in the case S{f, t) = S{r — vt)) 
and the ALH. And though we cannot at present provide general explanation of what 
makes such apparently different models be so closely connected, we hope that the results 
presented in Sections 01 El are rather convincing arguments in favour of the fact that this 
relation is useful, at least from the practical standpoint, as a tool for generating of a large 
number of solutions. On the other hand, this work presents 2D stationary structures of 
the isotropic continuous Heisenberg ferromagnetic spin system which have not been, to 
our knowledge, discussed in the literature and which seem to be interesting for the physics 
of magnetic phenomena. It should be noted that we have obtained our results in the 
framework of the classical model, and one of the most important questions that should be 
solved now, from the viewpoint of applications to magnetism, is to develop quantum, or 
at least semi-classical, theory of such structures. Another question we want to mention 
here is the following one. It is a widely known fact that solitons appear as a result of joint 
action of nonlinearity and some other mechanisms, such as dispersion. In our consideration 
we have neglected the temporal dispersion (temporal modulation), and its role has been 
played by the spatial one. So, it is interesting to take into account both temporal and 
spatial dispersions, because the competition of different mechanisms in nonlinear regime 
can lead to nontrivial results. These and some other related questions may be the subject 
of further investigations. 
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Appendix 



The aim of this section is to derive the matrix ^ related to the solution of equation (|1.1H1 
Shy 

s = m~^a^m (A.i) 

such that the matrices V ^ W 

y = 9^.^^i, W = d^-^-^ (A.2) 

have the structure of the ALH matrices 1)3.2(1 , (|3.3|) . The diagonalization ()A.1|) of a given 
matrix S is not unique. Suppose we have found a matrix F satisfying 

S = F-^a^F. (A.3) 
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Then any matrix 

^ = DF (A.4) 

with an arbitrary diagonal matrix D = diag {Dn, D22) satisfies ()A.1|) . Hence, to solve our 
problem we can start from any solution of ()A.3|) . for example from one given by 

F = l + a^S (A.5) 

(one can verify by simple calculations that this is indeed a solution of ()A.3|) ) and then to 
construct the diagonal matrix D such that the matrices 

V = d^ -^-^ = dD- D-^ + DdF ■ F'^ D'^, (A.6) 

W = 8"^ -^-^ = dD- D-^ + DBF- F-'^D-^ (A.7) 

possess the properties we need. 

To simplify the following formulae let us introduce the designation a, b and c for the 
elements of the matrix S, 

The main equation of this paper ()A.3|) can be rewritten now as a system 

a^s + aX + ^ {be' - b'c) = 0, (A.9) 

b^s + bX + i{ab' -a'b) =0, (A. 10) 

Czz + cX + i{ca' -c'a) =0, (A. 11) 



where 



X = a^a-z + -bzc-z + -bzC^, (A.12) 

/' = \^fz + \~^h (A.13) 



and 

,2 



a^ + bc=l. (A.14) 

Consider now the intermediate matrices V , W given by 



,-1 _ I I'll Vn 



Matrices V HA.6|) and W HA.7|) . which are the matrices V{X), W{X) from Section |21 can 
be written as 

^ ^^ ^^ ^ ' ^^-1/21 a In 1)22 + ^22 



D 



11 
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and 

\ L>ii 
Thus, if one takes the matrix D such that its elements satisfy 

(91nZ?ii = -Vii, 51nZ?22 = -V22, 

then matrices ()A.17|) . HA.18|) shah have the fohowing structure: 

y=fo *V w=(* ; 

As to other diagonal elements, they satisfy the identities 

5^22 = 9V22 -dW22, 
dWu =dWn-dVu. 
Using HA.9|) - (|A.1H) one can obtain that 

5^22 = dWii = — {b^Cz - bzCz) . 
4a 

From the definition of the matrices V, W, V, W and S it follows that 



^12^21 = V12V21 = -7 {al + 6.C,) 



,r2 
4 

and 



A.18) 

A.19) 
A.20) 

A.21) 

A.22) 

A.23) 

A.24) 
A.25) 

A.26) 
A.27) 

A.28) 
A.29) 

A.30) 

A.31) 
\ iXqo / 

i.e. to present them in the form 1)3. 2() . 1)3. 3|) . 

To summarize, we have derived, starting from a solution of the field equation (jl.llj) . 
the matrix ^, defined by HA.4|) . HA.5|) and HA.19|) . which can be used to perform the gauge 
transform from the Landau-Lifshitz linear problems to the ones of the ALH. 



W12W21 = W12W21 = -- (a? + bsc) . 
Using again ()A.9|) - I|A.11|) one can get that 

5yi2V^2i = -iX-'^dV22, 

dWi2W2i = -iX^dWu, 
i.e., comparing this result with HA.23|) . one can conclude that 

V'22 = -iX'^Vi2V2i+v{z), 

Wn = iX-^Wi2W2i + w{z). 
Thus, setting v{z) = —iX~^, w{z) = iX^ and defying the quantities qo, tq, qi, ri by 

qo = -^"^1^21, tq = iXVi2, qi = iXV2i, n = -iX'^Wu 

one can rewrite matrices (|A.17|) . ()A.18|) as 

/ -iA"Vo \ Y _ ( ^'^'^ ~ ^^0^1 iXri \ 

\ -iX'^qi -iA~2 + iroqi j ' I iXqo / ' 
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